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In a functional approach to QCD the infrared behaviour of Landau gauge Green functions is inves- 
tigated. Positivity violation for, and thus confinement of, gluons is demonstrated, and the analytic 
structure of the gluon propagator is studied. Quark confinement is related to an infrared diver- 
gence of the quark-gluon vertex. In the latter various components are dynamically induced due 
to the spontaneous breaking of chiral symmetry. As a result an infrared finite running coupling in 
the Yang-Mills sector is derived whereas the running coupling related to the quark-gluon vertex is 
infrared divergent. Based on a truncation for the quark-gluon vertex Dyson-Schwinger equation, 
which respects the analytically determined infrared behavior, numerical results for the coupled 
system of the quark propagator and vertex Dyson-Schwinger equations are presented. The result- 
ing quark mass function as well as the vertex function show only a very weak dependence on the 
current quark mass in the deep infrared. From this we infer by an analysis of the quark-quark 
scattering kernel a linearly rising quark potential with an almost mass independent string tension 
in the case of broken chiral symmetry. 
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1. Infrared behaviour of Landau gauge Yang-Mills theory 

Confinement, the Ua (1) anomaly, and dynamical chiral symmetry breaking are supposed to be 
properties of infrared QCD. Despite the progress achieved in our understanding of these phenomena 
the underlying mechanisms as well as possible interrelations between them are not yet uncovered. 
This talk aims at summarizing what we can learn from QCD Green functions in the Landau gauge 
about these topics. 

Explaining quark confinement, and hereby especially extracting the linearly rising static quark- 
antiquark potential and relating its properties to QCD degrees of freedom, has been the main ob- 
jective of many quite different studies. In ref. ^ some of these pictures where quark confinement 
is related to 

• the condensation of chromomagnetic monopoles (e.g. [^, ||]), 

• the percolation of center vortices(e.g. [|J]), 

• the AdSs / QCD correspondence (e.g. [||]), 

• the Gribov-Zwanziger scenario in Coulomb gauge (e.g. []6[ [7|]), or 

• the infrared behaviour of Landau gauge Green functions [|| ^, [I"0[ ], 

have been briefly reviewed. These explanations for confinement are seemingly different but there 
are surprising relations between them which are not yet understood. Given the current status one 
has to note that these pictures are definitely not mutually exclusive but simply reveal only different 
aspects of the confinement phenomenon. A similar statement on the relation between seemingly 
different scenarios is definitely true when considering the Ua(1) anomaly, see e.g. the discussion 



in ref. [11]. This should be kept in mind when in the following the results based on investigations 
of Landau gauge Green functions are presented. 

Employing Green functions dynamical chiral symmetry breaking reflects itself directly in 



terms of the quark propagator, see refs. M, RL 12, 13] and references therein. However, as we 



have demonstrated recently this picture is incomplete as it does not take into account the effects of 



dynamical chiral symmetry breaking on higher n -point functions []14[ |15|, |16p. This issue will be 
discussed in detail in this talk. To derive it one needs first some insight into the infrared behaviour 
of gluons and ghosts. 

1.1 Infrared exponents of gluons and ghosts 



q-l 




/q 

Figure 1: The Dyson-Schwinger equation for the ghost-gluon vertex. 

The possibly best starting point for the investigation of the infrared behaviour of gluons and 
jhosts is given by the Dyson-Schwinger equation for the ghost-gluon vertex function, see fig. [[]. 
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(NB: A MATHEMATICA package to derive Dyson-Schwinger equations can be found in [|17|].) In 
the Landau gauge the gluon propagator is transverse, and therefore one can employ the relation 



lfiD^ v (l-q) 



q^iD^ v (l-q). 



(1.1) 



to conclude that the ghost-gluon vertex stays finite when the outgoing ghost momentum vanishes, 
i.e. when — > [|18|]. This argument is valid to all orders in perturbation theory, a truely non- 
perturbative justification of the infrared finiteness of this vertex has been given in refs. Jl9| , [2C| , 

Knowing this property of the ghost-gluon vertex the Dyson-Schwinger equation for the ghost 
propagator, see fig. ^[ can be analysed. The only unknown functions in the deep infrared are the 




Figure 2: The Dyson-Schwinger equation for the ghost propagator, 
gluon and the ghost propagators, parameterized in the Euclidean domain as 



Z(k 2 
k 2 
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(1.2) 



In Landau gauge these propagators are best described by two invariant functions, Z(k 2 ) and G{k 2 ), 
respectively. As solutions of renormalized equations, these functions depend also on the renormal- 
ization scale ju. Furthermore, assuming that QCD Green functions can be expanded in asymptotic 
series, the integral in the ghost Dyson-Schwinger equation can be split up in three pieces: an 
infrared integral, an ultraviolet integral, and an expression for the ghost wave function renormal- 
ization. Hereby it is the resulting equation for the latter quantity which allows one to extract definite 
information [p2|] without using any truncation or ansatz. 



Recently it became clear that there are two distinct types of solutions, see ref. [23] and ref- 



erences therein. In this talk we will focus on the so-called scaling solution where the infrared 
behaviour of the gluon and ghost propagators is given by power laws. The corresponding expo- 
nents are uniquely related such that the gluon exponent is minus two times the ghost exponent 
As we will see later on this implies an infrared fixed point for the corresponding running coupling. 
The signs of the exponents are such that the gluon propagator is infrared suppressed as compared 
to the one for a free particle, the ghost propagator is infrared enhanced. 





+ 1/6 





+ 1/2 



Figure 3: The Dyson-Schwinger equation for the 3-gluon vertex. 
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Using this infrared power laws for the Yang-Mills propagators one can infer the infrared be- 
haviour of higher «-point functions. To this end the n-point Dyson-Schwinger equations have been 
studied in a skeleton expansion, i.e. a loop expansion using dressed propagators and vertices. Fur- 
thermore, an asymptotic expansion has been applied to all primitively divergent Green functions 



[25]. As an example consider the Dyson-Schwinger equation for the 3-gluon vertex which is dia- 



grammatically represented in fig. |[ Its skeleton expansion, see fig. f|, can be constructed via the 




•{JOT* + 

A Y » 5 T 



(a) (b) ^ (c) 1 (d) 

Figure 4: An example for the skeleton expansion of the 3-gluon vertex, 
insertions given in fig. H. These insertions have vanishing infrared anomalous dimensions which 

A T 

Figure 5: Insertions to reconstruct higher orders in the skeleton expansion. 

implies that the resulting higher order terms feature the same infrared scaling. Based on this the 
following general infrared behaviour for one-particle irreducible Green functions with 2n external 




ghost legs and m external gluon legs can be derived [25, 26]: 



T n ' m {p 2 ) ~ (p 2 )("- m ) K +( l -")( d / 2 - 2 ) (1_3) 

where K is one yet undetermined parameter, and d is the space-time dimension. 

Here two remarks are in order: First, exploiting Dyson-Schwinger and Exact Renormalization 



Group Equations one can show that this infrared solution is unique Q27Q. Second, there are addi- 
tional divergences when only some of the momenta of the n-point functions are vanishing [28]. 

1.2 Infrared fixed point of the running coupling in the Yang-Mills sector 

The infrared solutions ( |l.3[ ) include 

G( P 2 ) ~ (p 2 r K , z{ P 2 )~(p 2 Y K r 3 V)~(/r 3K , r*V) ~ (p 2 r 4K ■ (1-4) 

This implies that the running couplings related to these vertex functions possess an infrared fixed 
point: 

a^V) = «,GV)Z(p 2 ) - a 3 V) = a, [F^ (p 2 )] 2 z\ P 2 ) ~ S^L, 

a^(p 2 ) = a^(p 2 )Z 2 ( P 2 )~ C -^. (1.5) 
The infrared value of the coupling related to the ghost-gluon vertex can be computed to be [19|, 29]: 

a sk- g i m _ 47T r(3-2g)r(3 + *)r(i + K) 

[) ~6N C T*{2-k)T(2k) ' {b) 

This yields a**-* l (0) = 2.972 for N c = 3 and K = (93 - yi20T)/98 ~ 0.595353, which is the 
value obtained with a bare ghost-gluon vertex. 
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1.3 Gluon confinement by positivity violation for the gluon propagator 

Positivity violation for a propagator entails that the corresponding field is not related to an 
asymptotic state and thus to a particle. In the cases of gluons one can infer gluon confinement 
from such a scenario, see e.g. [ |30| ] and references therein. Therefore positivity violation of the 
propagator of transverse gluons has been for a long time a conjecture which has been convincingly 



verified by now, [ 31 , |32| ] as well as references therein. The basic feature is hereby the infrared 
suppression of transverse gluons caused by the infrared enhancement of ghosts. Being related to 
the confinement of tranverse gluons [ p0| ] it is certainly worth to have a closer look at the underlying 
analytic structure of the gluon propagator. 

As the infrared exponent K in the infrared power laws is an irrational number this implies 
already that the gluon propagator possesses a cut on the negative real p axis. It is possible to 
fit the solution for the gluon propagator accurately without introducing further singularities in the 



complex p plane [ pl| ] : 



Z &t (p 2 )=w( p2 ) (aUp 2 )V r . (1-7) 

Hereby w is a normalization parameter, and y = (—\3N c +4Nf)/(22N c — 4Nf) is the one-loop 
value for the anomalous dimension of the gluon propagator. The running coupling is expressed as 



[34]: 



. 2 . a s 0) An p 2 / 1 1 \ „ ON 

«fit P 2 = 7—^772 + ~R M — 2 W 2/A2 \ 2~7a2 r (L8) 

1 + p 2 / a qcd Po Aqcd + P 2 V Hr/Hco) P / a qcd - 1 / 



with j8o = (HA^c — 2A/y)/3. Note that the gluon propagator (1.7) possesses an important property: 
Wick rotation is directly possible! 

1.4 Summary on Yang-Mills sector 

To summarize this section we note that the here discussed scaling solution for Green functions 
in an Yang-Mills theory implies that the gluon propagator vanishes on the light cone, and n-point 
gluon vertex functions diverge on the light cone. Therefore an attempt to kick a gluon free, i.e. to 
produce a real gluon, immediately results in production of infinitely many virtual soft gluons to pro- 
duce perfect color charge screening. Positivity violation (which implies BRST quartet cancelation 



[J30t]) guarantees that this screening is total, or phrased otherwise, that gluons are confined. 



2. Dynamically induced scalar quark confinement 

2.1 Self-consistency between the quark propagator and the quark-gluon vertex 

As detailed above Landau gauge Green functions provide a consistent picture for gluon con- 
finement. However, due to the infrared suppression of the gluon propagator quark confinement 
seems even more unexplainable. To proceed it turns out to be necessary to study the Dyson- 
Schwinger equation for the quark propagator together with the one for the quark-gluon vertex 



, 10, 16, B3, 34]. Therefore a detailed study of this three-point function, and especially its 
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Figure 6: The Dyson-Schwinger equation for the quark-gluon vertex. 







Figure 7: The skeleton expansion for the quark-gluon vertex. 



infrared behaviour, is mandatory. Its Dyson-Schwinger equation is diagrammatically depicted in 
fig. ||, its skeleton expansion in fig. ^[ Hereby a drastic difference of the quarks as compared to 
Yang-Mills fields has to be taken into account: They possess a mass, and in addition one has to 
take into account that dynamical chiral symmetry breaking and thus dynamical mass generation 
may (and will) occur. 

The infrared analysis of the Yang-Mills theory has been generalized to full, although quenched, 



QCD in ref. [16]. Let us first consider a limit where the masses of the valence quarks are large, 
i.e. m > Aqcd- The remaining scales below Aq C d are those of the external momenta of the propa- 
gators and vertex functions. The Dyson-Schwinger equations are then employed to determine the 
selfconsistent solutions in terms of powers of the small external momentum scale p 2 -C Aqcd- The 
equations which have to be considered in addition to the ones of Yang-Mills theory are the ones for 
the quark propagator and the quark-gluon vertex which in turn have to be solved self-consistently. 

2.2 The quark-gluon vertex, chiral symmetry, and quark confinement 

The fully renormalized quark-gluon vertex consists of up to twelve linearly independent 
Dirac tensors. Half of these vanish in case chiral symmetry would be realized in the Wigner-Weyl 
mode, i.e. these tensor structures can only be non- vanishing either if chiral symmetry is explicitely 
broken by current masses and/or chiral symmetry is realized in Nambu-Goldstone mode (i.e. spon- 
taneously broken). From a solution of the Dyson-Schwinger equations we infer that these "scalar" 
structures are, in the chiral limit, generated non-perturbatively together with the dynamical quark 
mass function in a self-consistent fashion. Thus dynamical chiral symmetry breaking reflects itself 
not only in the propagator but also in the quark-gluon vertex. 

Performing an infrared analysis one obtains an infrared divergent solution for the quark-gluon 
vertex such that Dirac vector and "scalar" components of this vertex are infrared divergent with 
exponent — when either all momenta or when only the gluon momenta vanish [14, [TJ]. A nu- 
merical solution of a truncated set of Dyson-Schwinger equations confirms this infrared behavior, 
see fig. H. The driving pieces of this solution are the scalar Dirac amplitudes of the quark-gluon ver- 
tex and the scalar part of the quark propagator. Both pieces are only present when chiral symmetry 
is broken, either explicitely or dynamically. As can be seen from fig. ^| the function X\ multiplying 
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the tree-level tensor structure is only the leading one in the ultraviolet, in the infrared the scalar 
component A2 is even larger than the vector component X\ . 

b-quark 




-4 -3 -2 -1 1 2 

10 10 10 10 10 10 10 

P [GeV] 

Figure 8: The leading vector component Ai and a scalar component A2 of the quark-gluon vertex for a 
current mass corresponding to a b quark. Adapted from ref. 

The quark-gluon vertex may also serve to extract a running coupling. Using 

r*V) ~ iP 2 r 1/2 ' K , Zf{p 2 )~ const, Zip 2 ) ~ (p 2 ) 2K (2-1) 

one obtains 

a*V) = a, [V*(p 2 )] 2 [Z f ( P 2 )] 2 Z( P 2 ) ~ ^4 (2.2) 

and thus a coupling which is singular in the infrared contrary to the couplings from the Yang-Mills 
vertices. 

To determine whether this already relates to quark confinement the anomalous infrared expo- 
nent of the four-quark function is determined. The static quark potential can be obtained from this 
four-quark one-particle irreducible Green function, which, including the canonical dimensions, be- 
haves like (/? 2 )~ 2 for p 2 — > due to the singularity of the quark-gluon vertex for vanishing gluon 
momentum. Using a well-known relation for a function F <=c (p 2 )~ 2 yields 

V(r) = f-^F(p = 0,p)e^ ~ |r| (2.3) 

for the static quark-antiquark potential V(r). At this point one notes that, given the infrared diver- 
gence of the quark-gluon vertex as found in the scaling solution of the coupled system of Dyson- 
Schwinger equations, the vertex overcompensates the infared suppression of the gluon propagator 
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such that one obtains a linearly rising potential. The surprising fact is that this potential is dynami- 
cally induced and in the infrared dominantly scalar. 

To provide further understanding for the here found relation between chiral symmetry breaking 
and quark confinement one may keep chiral symmetry artificially in Wigner-Weyl mode, i.e. in the 
chiral limit one forces the quark mass term as well as the "scalar" terms in the quark-gluon vertex 
to be zero. The result of such a procedure is that then the running coupling from the quark-gluon 
vertex is no longer diverging but goes to a fixed point in the infrared similar to the couplings from 
the Yang-Mills vertices. Correspondingly, one obtains a l/r behaviour of the static quark potential. 
The "enforced" restoration of chiral symmetry is therefore directly linked with the disappearance of 
quark confinement. The infared properties of the quark-gluon vertex in the "unforced" solution thus 
constitute a novel mechanism that directly links chiral symmetry breaking with quark confinement. 

2.3 Mass dependence 

As mentioned above, based on a truncation for the quark-gluon vertex Dyson-Schwinger equa- 
tion which respects the analytically determined infrared behavior, numerical results for the coupled 
system of the quark propagator and vertex Dyson-Schwinger equations have been obtained []i6|]. 
In fig. ^ the resulting quark mass function is presented for the values of current masses related to 
u/d,s,c and b quarks. These quark mass functions as well as the vertex functions show only a very 



icVe — — — — d 




2 



P [GeV] 

Figure 9: The quark mass functions for several values of current masses. Adapted from ref. [ |l6| | 

weak dependence on the current quark mass in the deep infrared. First of all, the infrared exponent 
for all vertex functions are mass independent. Second, comparing the coefficients of these power 
laws, given in table |T], we see a much smaller dependence on the quark mass than could have been 
anticipated from a naive analysis. As the static potential will be build up by the infrared divergent 
parts of the quark-gluon vertex we thus predict that the string tension is almost mass-independent. 
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u/d 


s 


c 


b 


m(n 2 )[MeV] 


3 


50 


1200 


4200 


M{0)[MeV] 


270 


270 


320 


650 


X^GeV 1 ! 2 ^} 


3.95 


3.60 


2.00 


2.73 


X 2 [GeV l / 2+K ] 


8.70 


7.97 


4.44 


6.08 



Table 1: The infrared coefficients for the dressing functions A12 for different current quark masses m(jx 2 ) (at 
the renormalization point fx 2 — 170GeV). Also given are the values of the mass function at zero momentum. 

Provided the Coulombic part of potential has still a large effect in the bottomonium spectrum this 
is in agreement with experimental splittings of charmonia and bottonia systems. 

Thus our analysis provides (at least) two surprises: First, the longstanding discussion on 
whether confinement is of a vector or of a scalar nature is oversimplified. Dynamically gener- 
ated structures of vector, scalar and tensor types in the quark-gluon vertex lead to a rich structure of 
the quark-antiquark interaction. Second, the mass dependence of the string tension does not seem 



to follow any easy rule, cf. also recent improved lattice data [p5|]. These and other details of the 
confining interaction clearly deserve more investigations. 

2.4 Summary and outlook on quark sector 

The scaling solution for Landau-gauge QCD Green functions has far reaching consequences 
for the infrared behaviour of quarks. The resulting dynamical chiral symmetry breaking and the 
related dynamical mass generation lead to an infrared trivial quark propagator. But dynamical 
chiral symmetry breaking also occurs in the quark-gluon vertex. Hereby several components, and 
especially also the "scalar" ones, diverge on the quark "mass" shell. 

An attempt to kick a quark free, i.e. to produce a real quark, immediately results in production 
of infinitely many virtual soft gluons. Hereby these gluons do not only couple vector-like but also 
scalar-like! And the vertex function diverges such that effectively a linearly rising potential is 
produced. One obtains therefore infrared slavery and quark confinement. 

Although these results are in itself quite surprising and encouraging they provide only the start- 
ing point for an understanding of quark confinement in functional approaches. Questions which 
will be investigated further include: 

• How does the formation of a string between colour sources reflect itself in Green functions? 

• Can we learn about the properties of the confining field configurations from Green functions? 

• How does all this relate to other (partially) successful approaches to understand the phe- 
nomenon of quark confinement? 
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